Z2 topological zero-energy modes in commensurate Aubry- Andre-Harper models. 



Sriram Ganeshan,^ Kai Sun,^'^ and S. Das Sarma^ 

^Condensed Matter Theory Center and Joint Quantum Institute, 
Department of Physics, University of Maryland, College Park, MD 20742, USA 
^Department of Physics, University of Michigan, Ann Arbor, MI 48109, USA 

(Dated: January 25, 2013) 

Aubry-Andre (AA) model has been the subject of extensive theoretical research in the context 
of quantum localization. Recently, it is shown that one-dimensional quasicrystals described by the 
incommensurate Aubry-Andre model has non-trivial topology In this paper, we show that the 
commensurate off-diagonal Aubry-Andre model is topologically nontrivial in the gapless regime and 
supports zero-energy edge modes with Z2 index. Unlike the incommensurate case, the nontrivial 
topology in the off-diagonal Aubry-Andre model is attributed to the topological properties of the 
one-dimensional Majorana chain. We discuss the feasibility of experimental observability of our 
predicted Z2 topological phase. 



Introduction. Anderson localization is a quantum- 
interference-induced disorder-tuned quantum phase tran- 
sition on a tight binding lattice where the system wave- 
function changes from being extended ("metal") to ex- 
ponentially localized ("insulator") at a critical value of 
the disorder strength [Ij. In one-dimensional (ID) sys- 
tems, Anderson localization is trivial since the critical 
disorder is zero, and all states for any finite disorder are 
localized. The absence of a true quantum phase transi- 
tion makes ID Anderson localization rather uninterest- 
ing from the perspective of the physics of disorder-tuned 
metal-insulator transition. However, Aubry and Andre 
(AA) [2 predicted the existence of a localization tran- 
sition for certain ID quasiperiodic systems, where the 
transition arises from the existence of an incommensu- 
rate potential of finite strength mimicking disorder in a 
ID tight binding model. The form of this quasiperiodic 
potential is usually chosen to be a cosine function which 
is incommensurate with the underlying periodic tight- 
binding ID lattice. This result has led to extensive theo- 
retical studies of disorder effects in the AA model during 
the last few decades [3]-[7]- Recent experiments [H |9] 
have realized the quasiperiodic AA model in optical lat- 
tices and have observed the signature of a localization 
transition [8 in agreement with the theory [2 . 

One very interesting aspect of the ID AA model is 
that it can be exactly mapped to the 2D Hofstadter 
model [To', TT. The Hosftadter model describes the 
topologically-nontrivial 2D quantum Hall (QH) system 
on a lattice |T2UT5] . This mapping implies that the ID 
AA model must have topologically-protected edge states 
similar to the gapless edge states of the QH effect. Re- 
cently, these edge states have been observed experimen- 
tally [9] and this mapping has been used to topologically 
classify ID quasicrystals described by the incommensu- 
rate AA model [H [TBHIS] . In this letter, we study the 
commensurate off-diagonal AA model. We note in the 
passing that although we use the AA nomenclature to 
discuss our model (mainly to establish connection with 
existing work in the literature), our proposed commen- 



surate off-diagonal system is simply a 'bichromatic' ID 
system with two competing underlying ID commensu- 
rate periodic potentials with arbitrary phases. The hop- 
ping amplitude of the off-diagonal AA model has a co- 
sine modulation in the real space commensurate with the 
lattice. This is in contrast to the diagonal AA model 
which has a cosine modulation in the potential energy 
term. Both these versions can be unified within a gen- 
eralized A A model (also known as generalized Harper 
model [ini HSl [19] - - we use the AA model terminol- 
ogy throughout this paper). In particular, we focus on 
the parameter range where the AA model is gapless and 
thus cannot be mapped onto a QH system. Surprisingly 
we find that edge states exist in this seemingly topo- 
logically trivial model. These edge states are topologi- 
cally protected, but they belong to a different topolog- 
ical class than the QH edge. The topological origin of 
these edge states are analogous to the zero-energy edge 
states along the zigzag edge of graphene [2QH22] and are 
directly connected to the Z2 topological index of the Ki- 
taev model [23]. To make the Z2 topology transparent, 
we rewrite our model in the Majorana basis. It must 
be emphasized that the commensurate off-diagonal AA 
model studied in this work has extended Bloch ID band 
bulk states for all parameter values (and no localization 
transition at all), and our establishing its topological Z2 
edge behavior differs qualitatively from all earlier recent 
work on topological properties of the incommensurate 
AA model ^ [MMl [21 [25]. 

In addition, we provide a specific experimental setup 
for the realization of the off-diagonal AA model. We 
propose a ID lattice composed of coupled single mode 
waveguides with varying lattice spacings mimicking the 
commensurate modulation of the hopping strength. To 
understand whether the topological edge states can be 
observed in a real experimental setup, we examine the 
robustness of these edge states against next order hop- 
ping and fluctuations in lattice potentials, both of which 
arise in a real experimental system. Although the Z2 
index cannot be clearly defined in the presence of these 
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(a)Topologically-nontrivial off-diagonal AA model for (b)Topologically trivial diagonal AA model for 6 = ^ 

6=1 (tt-Aux) (tt-Aux) 



Figure 1. Upper panel shows the two zero-energy eigenstates (red and blue lines) plotted as a function of position for three 
different values of The lower panel is the energy spectrum plotted as a function of Lp for 100 sites, (a) Spectrum and 
eigenstates with parameters A = 0.4, v — ^ and h — 1/2. (b) Spectrum and eigenstates with parameters v — 0.4, A = and 
h= 1/2. 



realistic effects, explicit numerical calculations indicate 
that these edge modes are robust and observable for a 
wide range of parameter space. This robustness, arising 
from the topological nature of the zero modes, is of cru- 
cial importance for the experimental verification of our 
prediction. 

Model. We consider the generalized ID AA model, 
which is described by the following Hamiltonian 

N-l 

H = ^ t[l + Acos(27r6n + (px)]cl^-^Cn + h.c. 

N 

+ ^ V cos(27r6n + (py)clcn. (1) 

n=l 

This ID chain has N sites (n = 1, 2, . . ., N). We adopt 
open boundary conditions with n = 1 and n = N be- 
ing the two edge sites. To be consistent with previous 
literature on topological edge modes [23 , we consider 
fermionic particles which are created and annihilated by 
fermioinic operators cj^ and c^. We emphasize, however, 
that our work and all conclusions are equally valid for 
the corresponding bosonic case since we are considering 
a noninteracting ID quantum system. The first term in 
the Hamiltonian is the kinetic energy from the nearest- 
neighbor hopping, and the last term describes the on-site 
potential energy. The inhomogeneity in hopping strength 
and potential energy terms are described via cosine mod- 
ulations of the strength A and v respectively. The cosine 
modulations have periodicity 1/b and phase factor cpx 



and (fy. The special case with X = {v = 0) corresponds 
to the diagonal (off-diagonal) AA model. The general- 
ized AA model can be derived starting from an ancestor 
2D Hofstadter model with next-nearest-neighbor (diag- 
onal) hopping terms [16l [191 [2S]- Starting from a 2D 
ancestor, the phase terms (fx and (fy correspond to the 
same higher dimensional momentum degree of freedom 
and are related by (fx = -\- ^rb. Experimentally, one 
can design setups where both cpx and (fy can be tuned 
independently. To make contact with the experimental 
results we keep our notations general with cpx and cpy as 
independent variables. 

For irrational 6, the diagonal AA model (A = 0, 7^ 0) 
shows a localization transition as v is increased beyond 
the critical value {v = 2t) with all states being ex- 
tended(localized) for v < 2t {v > 2t) For 
rational 6, it is known that by treating the phase (fy as the 
momentum of another space dimension, the diagonal AA 
model can be mapped onto a 2D Hofstadter lattice with 
27r6 magnetic flux per plaquette [TOl [H]. For b ^ 1/2, 
the Hofstadter lattice has gapped energy bands with non- 
trivial topology, described by nonzero Chern numbers. 
Due to the nontrivial topology, localized edge modes are 
expected for a finite-sized system with boundary. It is 
worthwhile to mention here that although the mapping 
to Hofstadter lattice is well-defined for rational values of 
6, the topologically-protected edge states remain stable 
even if b takes irrational values. We emphasize here that 
the topological equivalence between the diagonal and off- 
diagonal AA model exists only for the irrational values 
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of b (incommensurate case) [9l [16] . 

For the special case of b = 1/2, which is the focus 
of this work, the AA model can be mapped onto a 2D 
Hofstadter model with tt flux per plaquette. Under time- 
reversal transformation, a tt flux simply turns into a — tt 
flux. Since the magnetic flux terms are only well defined 
modulo 27r for a lattice, the system is then invariant un- 
der the time-reversal transformation. Therefore, the sys- 
tem shows no QH effect and thus has no QH edge modes. 
In fact, this Hofstadter lattice model has no band gap 
but contains two Dirac points in analogy to graphene. 
This gapless 7r-flux state has been extensively studied 
in the context of algebraic spin liquids (see for example 
Refs. [28l|29] and references therein). 

We have calculated the energy spectrum for diagonal 
and off-diagonal AA models at different values of (fy and 
(fx (Fig.jl]). For periodic boundary conditions (not shown 
here), both models show two energy bands with two Dirac 
points. If we perform the same calculations on a ID lat- 
tice with open edges, the energy spectrum (as a function 
of (px^v) for the diagonal A A model remains the same as 
in the case of periodic boundary conditions. But for the 
off-diagonal AA model, zero-energy states are observed 
for — 7r/2 < (fx < 7r/2 as shown in Fig. [l] By examin- 
ing the wavefunctions associated with these zero-energy 
states, we find that they are actually boundary states lo- 
calized around the two edges of the system, as shown in 
the inset of Fig. [l] These are the topological zero-energy 
(edge) modes alluded to in the title of our paper. 

Degenerate Majorana modes. The edge states we find 
here have a topological origin, which can be understood 
analytically by rewriting the off-diagonal AA model (for 
b = 1/2) in the Majorana basis. We define C2n = 72n + 
ir2n, C2n+i = T2n+i + ^72n+i, whcrc 7 and r are two 
species of Majorana fermions. In this new basis, the off- 
diagonal A A model becomes, 

H = ^ A_(72n-l72n) + A+ (72n72n+l ) 
n 

- ^ A_(r2n-ir2n) - A+(r2nr2n+i), (2) 

n 

where A± = 2zt(l±2A cos cpx)- Here, the system contains 
two identical ID Majorana chains which are decoupled 
from each other. In the study of ID topological super- 
conductors, it is known that a Majorana fermion chain 
supports a Z2 topological index [23J. For |A+| > |A_|, 
the Majorana chain is topologically nontrivial and has 
one zero-energy Majorana mode localized at each edge. 
For the opposite regime | A+| < | A_ |, the system is topo- 
logically trivial with no edge modes. For our model, this 
implies that two Majorana modes, which are equivalent 
to a Dirac edge mode, are expected at each of the two 
edges for cos cpx > 0^ which agrees perfectly with the nu- 
merical results shown in Fig. [T] For odd number of sites, 
there exists a single Majorana mode localized on one of 
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Figure 2. Energy spectrum with odd number of sites (N — 
101). Single zero-energy mode is always localized on either 
one of the edges except for the Dirac points as shown by the 
wavefunction plots. 



the edge sites except for the Dirac points as shown in 
Fig. (§. 
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Figure 3. (a) Energy spectrum with next nearest neighbor 
hopping for N = 100 sites (t = 1, = 0.2, A = 0.4 and 
= 0). (b) Degeneracy of the zero modes is lifted with weak 
onsite modulation(t — 1^ v — 0.2 and A = 0.4 ) 



Robustness. We now test the robustness of these edge 
modes. As long as the particle-hole symmetry [cn 
( — l)^cj^ and cj^ (—1)^0^] is preserved [20 , the map- 
ping to two decoupled Majorana chains remains valid and 
thus the edge modes are stable with their energy pinned 
to zero. However, it is worthwhile to point out that 
such a symmetry can be broken explicitly by the next- 
nearest-neighbor hopping and a modulating potential en- 
ergy, which will couple the two Majorana fermion chains 
together. This coupling will result in the Z2 topological 
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index being ill- defined. To test the fate of the edge modes 
under such relevant perturbations, we introduce next- 
nearest-neighbor hopping and a modulating potential en- 
ergy to the off-diagonal AA model. For simplicity, we as- 
sume that the strength of the next-nearest-neighbor hop- 
ping is site independent {Hnnn = Zl^ + h.c). 
The modulating potential energy is introduced using the 
onsite term v in Eq. ([T]) and without any loss of gener- 
ality we assume cpy = (fx = cp from hereon. As shown 
in Fig. |3| the edge modes are found to be very robust. 
Although the energies are shifted away from zero, the lo- 
calized edge modes (see wavefunction plot in Fig. |3| can 
be observed even for very strong f and v. This robust- 
ness is a direct result of the topological nature of their 
origin. 
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Figure 4. Tight binding energy spectrum with open boundary 
conditions for 100 sites and b = l/4.(a)Energy spectrum as 
a function of the parameter cp with A = 1.0 and v = 0. (b) 
Central gap closes for A = V2. (c) Energy spectrum with 
A = 3.0. Position of degenerate zero modes shifts by 7r/2 (d) 
t' — 0.2 opens up a gap in the bulk with gap less edge modes 
of QH type connecting the bulk bands. 

Furthermore, the above analysis can be generalized for 
the case of 6 = 1/(2(7), where is a positive integer. For 
example, the off-diagonal AA model with 6 = 1/4 has 
four energy bands. The top and bottom bands are fully 
gapped (from the two central bands) but the two bands in 
the middle have four band crossing (Dirac) points located 
at (p — 71 /^^ 37r/4, 57r/4 and 77r/4, as shown in Fig. |4] 
Inside the band gap between the top two bands (or be- 
tween the bottom two bands), the quantum Hall edges 
states can be observed. Between the two central bands, 
zero-energy edge states are observed. For \/t < \/2, the 
zero-energy edge modes are found for J <(/:)< ^ and 
^ < if < ^ . When X/t > V^, the zero modes are found 
for — ^ < (f < J and ^ < (f < The marginal case 
with X/t = \/2 shows no gap at any value of Lp and thus 
has no edge states between these two bands. These zero- 
energy states are of the same origin as the zero-energy 



edge states discussed above for 6 = 1/2. Writing 6=1/4 
case in the Majorana basis we obtain 

a = ^ A2n-l(72n-l72n) + (72n72n+l ) 

n 

- ^ A2n-l(T2n-l'r2n) " ^2n{^2n^2n+l) , (3) 
n 

where = 2it[l + Xcos{(p + )]• The zero-energy 
modes exist for the parameter range satisfying |1 — 
A cos (p\ > |1 — A sin (p\ which is in prefect agreement with 
the numerical results. For q > 1 the next nearest neigh- 
bor hopping term {t') or diagonal cosine modulation {v) 
opens a gap between the two central bands, and the zero- 
energy edge modes adiabatically turn into mid-gap edge 
modes. (Fig. [4jd). 
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Figure 5. Experimental schematic: The red lines shows the 
single mode waveguides. The hopping amplitude is modulated 
by varying the spacing between the waveguides. The phase (px 
can be scanned along the length of the waveguide by smoothly 
changing the spacing between successive channels. 



Experimental realization. Recent experimental ad- 
vances [H E] have been successful in realizing topologi- 
cal features of ID AA models for incommensurate lattice 
modulations. Topological properties of the half flux state 
can be realized in photonic crystals using similar setups 
as used in Ref. [9 . As shown in Fig. [sj each waveg- 
uide corresponds to a site in the tight binding Hamil- 
tonian. The potential energy at each site can be ad- 
justed by tuning the width of each waveguide, which of- 
fers a way to control v and cpy. The off-diagonal mod- 
ulation A and cpx can be controlled by tuning the spac- 
ing between two neighboring waveguides, which deter- 
mines the hopping strengths. By varying these spac- 
ings, one can tune the value of cpx into the topological 
regime — 7r/2 < (/:?a < 7r/2 to observe the topological edge 
modes discussed above. In addition, one can scan dif- 
ferent value of cpx by smoothly varying the spacing be- 
tween the waveguides along the direction of the waveg- 
uides. In principle, it should also be possible to experi- 
mentally study our proposal by using suitably designed 
cold atom optical lattices [30l ISj or semiconductor struc- 
tures [32l [33] where ID AA-type quantum systems have 
earlier been realized in the laboratory. 
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Conclusion. In this work we have unearthed a novel 
topological aspect of the commensurate off-diagonal AA 
model (i.e. ID bicromatic lattice model). It is shown 
both analytically and numerically that the 6 = 1/2 flux 
state of the off-diagonal AA model supports a topo- 
logically non-trivial zero-energy edge mode that has a 
Z2 topological index with respect to the higher dimen- 
sional phase parameter ipx. The Z2 topological nature of 
the commensurate bichromatic ID Aubry- Andre-Harper 
model uncovered by us shows some deep and surprising 
connections between simple ID hopping models on the 
one hand with spin liquids, graphene, and ID topologi- 
cal superconductivity on the other. 
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